Consensus time and conformity in the adaptive voter model 

Tim Rogers 
Department of Mathematical Sciences, Bath University, Claverton Down, Bath, BA2 TAB, UK 



cn 

o 

(N 

< 



Oh 
o ■ 

o- 
o 

• I— I 

>.. 

^• 



> 

O 
m 



X 




Thilo Gross 
Department of Engineering Mathematics, University of Bristol, Woodland Road, Bristol, BS8 1 UB, UK 

The adaptive voter model is a paradigmatic model in tiie study of opinion formation. Here we 
propose an extension for this model, in which conflicts are resolved by obtaining another opinion, 
and analytically study the time required for consensus to emerge. Our results shed light on the rich 
phenomenology of both the original and extended adaptive voter models, including a dynamical 
phase transition in the scaling behavior of the mean time to consensus. 

PACS numbers: 89.75.Fb 



In nature, collective intelligence is observed in a wide 
variety of species. Quite generally groups of animals are 
able to aggregate information and make decisions col- 
lectively [l|. The most impressive example is perhaps 
human culture, which is created through the aggrega- 
tion and transmission of individual insights and opinions. 
However, while collective decision making seems to be 
universally beneficial in animals, it can have adverse ef- 
fect in humans, where the exchange of opinions can lead 
to the propagation of counter-factual rumors and beliefs, 
fuels youth violence, and can even give rise to the for- 
mation of radicalized groups. We thus have to ask under 
which conditions collective information processing leads 
to the discovery of new factual truth and when it results 
in counter-factual myths. In particular in the light of the 
profound changes that the opinion formation process in 
humans is currently undergoing, a deeper understanding 
of this process is desirable, to exploit the beneficial effects 
of new technologies and mitigate detrimental ones. 

A paradigmatic model in this field is the adaptive voter 
model [2|, [Sj , describing a collection of individual agents 
whose opinions and social contacts may change over time. 
Agents hold one of two opinions, say A and B, and are 
linked together by a sparse network of social interactions. 
The system evolves in time as follows: pairs of connected 
agents with opposing opinions are randomly chosen and 
either (i) the conflict is resolved by one agent adopting 
the opinion of the other, or (ii) one agent breaks the 
contact and forms a new link to a different agent. After 
a sufficiently long time the system reaches one of two 
types of absorbing state: a consensus state in which all 
agents hold the same opinion, or a fragmented state in 
which both opinions survive in disconnected components 

While the adaptive voter model has been explored in 
several recent studies, a larger body of previous work fo- 
cuses on opinion dynamics on static networks (where the 
rewiring process does not occur and hence fragmentation 
is impossible). The main question addressed in these 
studies is the time taken for consensus to emerge [7h9|. 
Generally, the formation of consensus is slow, with the 



FIG. 1. (Color Online) Dynamical phase transition in the 
scaling behavior of the mean time to reach consensus Tc as a 
function of network size N. At the critical value p = 1/2 the 
growth of Tc is linear, whilst being exponential for p above 
the critical point, and logarithmic below. In each case circles 
give the average over 100 simulation runs with fc = 10 and 
(j) = 0.1, while the solid lines show the result of the theory 
developed in the main text - see equation @. 



time taken growing like N^, where the exponent /i < 1 
depends on the degree distribution of the underlying net- 
work 9]. In contrast to this work, the major results in 
the adaptive networks literature have all concerned the 
occurrence of fragmentation [3h6| , whereas the question 
of consensus time has so far been largely neglected. 

Here we show analytically that when consensus is 
reached in the adaptive voter model, the time taken 
grows linearly with N. Furthermore, we show that when 
we allow pairs of agents to resolve conflicts by seeking 
another agent's opinion, this can either speed up or slow 
down the formation of consensus. This extension of the 
adaptive voter model exhibits a dynamical phase transi- 
tion between exponential, linear, and logarithmic scaling 
laws, depending upon the probability of accepting the 
other agent's opinion. 

Consider a network of N nodes (agents) joined by a 
total of K edges, which represent social interactions. Ini- 



tially each agent is randomly assigned either opinion A or 
B, and the locations of the edges are placed randomly. In 
each time step an edge {i,j) is selected at random. If the 
focal edge connects agents that hold different opinions 
then it said to be active and the corresponding conflict 
is resolved in either of two ways. With probability (p the 
edge is rewired, with node i cutting the edge and creating 
a new one to another node, selected at random from the 
set of all nodes holding the same opinion as i. If the edge 
is not rewired then a third opinion is sought: another 
node is selected at random from the rest of the graph 
and with probability p both i and j adopt the opinion of 
the new node, otherwise, both i and j adopt the opposite 
opinion. 

The parameter p can be interpreted as a measure of 
social conformity [lO| and may range from w 1 for a 
strongly conformist opinion formation process, to ~ 
for nonconformist processes. For p — 1/2 the decision is 
not biased by the other agent, and the standard adap- 
tive voter model is recovered. In most relevant contexts 
the stochastic response to another agent should not be 
interpreted as an actual consultation, but rather as an 
influence from the cultural 'mean field' propagated by 
mainstream mass media. 

As we will see, the model exhibits a dynamical phase 
transition as p crosses the threshold pc = 1/2. Simulation 
results (Fig. 1) show that when agents accept the third 
opinion with probability p > pc, the time to consensus 
only grows logarithmically with N, whereas in the case 
that p < Pc it exponentially diverges. At precisely p — Pc 
the original adaptive voter model is recovered, and we 
find a linear growth of the consensus time. 

Our analysis focuses on the time evolution of the first 
and second order moments of the distribution of node 
states. We denote the numbers of agents with each opin- 
ion by [A] and [B], and the numbers of edges between 
different agents by [AA], [AB] and [BB]. Because of 
the conservation laws for nodes, [A] + [B] = N, and 
edges, [AA] + [AB] + [BB] = K, the state of the system 
can be summarized by just three independent quantities: 
X =[A]-[B],Y = [AA] - [BB] and Z = [AB]. In anal- 
ogy with spin glasses, the first two of these describe the 
magnetization of nodes and edges, while the third speci- 
fies the number of active edges and therefore controls the 
overall reaction rate of the system. 

For understanding the role of demographic noise in 
opinion formation it is necessary to have a fully stochas- 
tic description of the time evolution of our order parame- 
ters. We achieve this by defining a Markov jump process 
which acts as a pair based proxy (PBP) for the adaptive 
network model [ll|. Introducing the system state vector 
il = {X, Y, Z)"^ , we enumerate the possible events which 
alter the system state. For each of these events, we com- 
pute the probability ri of occurrence in a given time step, 
and the average net change Si to il caused by the event. 
These calculations may be performed explicitly, however 



they require a little thought; full details are given in the 
supplementary material. 

The PBP model is then defined as follows: in each time 
step a jump vector Si is chosen randomly with probability 
ri , and the summary vector is updated by Jl i-> Jl + Sj . 
The PBP represents a considerable reduction in complex- 
ity from the original adaptive network, and yet retains 
the essential stochastic nature of the system. 

In the limit of large network size, the PBP can be 
further reduced to a low-dimensional system of SDEs. 
For simplicity, we package the update vectors into a sto- 
ichiometric matrix S — (si • • • S4) and collect the event 



probabilities in a vector r 



1^1 



■ri 



Defining the 



rescaled variables x = XJN , y — Y/ K and z/ K and we 
apply Kurtz' theorem 12] to obtain the following SDE 
in the limit of large N: 




Sr + —= rjit) . 



where 



{V^{t)flJit'))^Sit~t')B,J, 

and the noise correlation matrix B is given by 
Bij — 2_^ SikTkSkj ■ 



(1) 



(2) 



(3) 



When written explicitly in terms of x, y and z (see sup- 
plement), the equation for the magnetization is the most 
instructive: 



^=4(l-</))(p-l/2)a:z+^ryi(i) 
di Jn 



(4) 



The factor of {p — l/2)x constitutes either positive or 
negative feedback depending on the value of p, which 
suggests a transition in behavior around point pc — 1/2. 
The nonlinear interaction xz shows that the for the dy- 
namics to progress requires the presence of active edges, 
as well as an overall imbalance of opinions. 

Wc make analytical progress by exploring the behavior 
of ([T]) in the neighborhood of the transition, introducing 
£ — p — 1/2. Let us first consider the case e = in the 
deterministic limit N —> 00. In this limit, the system ([T]) 
possesses two manifolds of fixed points. The first is the 
plane z = 0, which represents the state in which there 
are no active edges and thus the dynamics are frozen, 
note that these are also absorbing states of the finite- 
size network model, corresponding to fragmentation. The 
second manifold of fixed points defines a parabola: 



y 



1(1 -x^)^ 
2^ ^ 1 



(5) 



with (j)^, = {k — 2)/(fc — 1) being the approximate critical 
rewiring rate for the fragmentation transition identified 



in y . We note that Eq. ([5]) is a pair- level approximation 
which becomes poor close to the fragmentation point; 
see Fig. [2b- However, the question of consensus time 
concerns only values of (j) below the transition, where 
we find the approximation to be sufficient for a large 
parameter range. 

Further insight is gained from a linear stability analy- 
sis. We write J for the Jacobian matrix of (jlj. At every 
point on the parabola ([5]) , we find J has eigenvalues 

Ai=0, A2-0-</.*, \,= ^^^-^*\ (6) 

Z — 0* 

The corresponding eigenvectors are 



Vl 




,V2 




,«3 




X 1 


(7) 


l/(fc-l)(l- 
The second 



where the constants are given hy fii — — 

(/)) and /i2 = -1 - 2/(fc - l)(fc - 2)(1 - 

two eigenvalues are negative, and large in comparison to 

Ai = 0, meaning that trajectories close to the parabola 

collapse quickly in the directions of V2 and 1)3. We show 

an example of this behavior in Fig. [^b; it is key to the 

remaining analysis. 

Although the full stochastic system ([Ij cannot be 
easily solved, we can derive an 'effective' solvable one- 
dimensional system by restricting our attention to behav- 
ior in the neighborhood of the slow manifold. We reason 
as follows: in short time windows small Gaussian pertur- 
bations described by the noise correlation matrix defined 
in ([3]) may move the system away from the slow manifold; 
for sufRciently small perturbations the net drift is then 
governed by the fast eigenvectors 1)2,3. A mathematical 
derivation of the procedure (in analogy to 13|, llJ]) is 
given in the supplementary material. 

We obtain a reduced equation for motion on the man- 
ifold, in which x is the only remaining variable. 



Ax 
di 



2e(0* - (j))x{l - x^) + 



2(0, -(/.)(l-a;2) 



N 



?w 



(8) 

where ^ is a standard Gaussian white noise variable. The 
picture we have now is as follows: from the initial condi- 
tion the system state moves rapidly to the parabolic slow 
manifold ([5]), where it then drifts stochastically accord- 
ing to ([U until eventually reaching one of the absorbing 
consensus states at x = ±1. 

We are interested in the mean waiting time before con- 
sensus is reached. Since our theory is one-dimensional, 
we follow 151 to derive 



Tc^ 



N 



y pSN{x''-y'') 



■ dxdy . 



(9) 



/O Jo 1 — u. 

In the special case e — Q the integral above can be com- 
puted easily to obtain 

log(2) 



Tn 



N- 



(10) 




FIG. 2. (Color Online) (a) Dependence of consensus time Tc 
on the rewiring rate cj>, when e = and k = 10. Orange 
circles show the average of 100 samples, the solid black line 
is the theoretical prediction (|10p . while the dashed line in- 
dicates the point of the fragmentation transition, as derived 
in [5|. (b) Comparison between the slow manifold ((5]) and 
typical simulation trajectories for (j) — 0.5 (red) and (j) — 0.7 
(blue). The discrepancy between simulations and theory in 
the case (p — 0.7 illustrates the breakdown of the pair-level 
approximation close to the fragmentation transition. 



In Fig. [2] we show a comparison between this prediction 
and the results of numerical simulations; the agreement 
is excellent for values of (f) away from the fragmentation 
transition. 

If p < 1/2 then e is negative and the large N asymp- 
totic of (O can be computed by Laplace's method, yield- 



ing 



,\e\N 



Tc 



4(0, -0)V \e\^N 



(11) 



Alternatively, for p > 1/2 the system ([8]) is determin- 
istically unstable and thus the main contribution to Tc 
comes from the initial symmetry breaking perturbation. 
In the neighborhood of the initial condition a; = 0, we 
have the linearized equation 



dx 
di 



2e{(j).^, - 4>)x - 



N 



■m- (12) 



Following [3 again we find {x^) = {e^<'l>'-'P)t - l)/2eiV, 
and thus the time taken for x to reach a magnitude of 
order one is 



IokA^ 



4e(0, - 0) 



(13) 



A comparison of these predictions with numerical results 
(Fig. [3|) shows an excellent agreement. These results es- 
tablish a trichotomy between exponential, linear and log- 
arithmic scaling laws, dependent on the parameter p. 

Let us now discuss the sensitivity of the model to small 
changes in the microscopic rules. In a recently proposed 
variant of the adaptive voter model the target nodes in 
rewiring events are chosen randomly without regard to 
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FIG. 3. (Color Online) Large A'' scaling for p either side of 
the critical value pc = 1/2. On both plots circles show the 
average consensus time for 100 simulation runs with A; = 10 
and (/> = 0.1. On the left the black line is given by Eq. (|ll|l . 
while on the right the slope is given by Eq. (|13p . whereas the 
intercept has been chosen manually for comparison. 



their opinion (e.g. [6|). We refer to this variant as rewire- 
to-random, as opposed to the rewire-to-same scheme we 
considered above. This change in model specification re- 
sults in a different expression for the event probabihty 
vector r, however, the rest of the analysis may be re- 
peated without difficulty (see supplement). The result- 
ing expressions are quite complex, however, expanding in 
large k we find the simple form 



Tc 



N 



log(2) 



0{k- 



(14) 



Comparing to (TTO)) we see that, for sufficiently large fc, 
the choice of rewiring rule does not much change the 
typical time to reach consensus. This prediction is tested 
numerically in Fig. |4^. 

The mechanism for choosing nodes to update may also 
be altered. We have so far considered link update, mean- 
ing that at each time step an edge is chosen at random 
and then one of the nodes connected to that edge is se- 
lected to update. Alternative model formulations use 
node update rules [la], where one first chooses a node i 
before selecting one of its neighbors j, and then in direct 
node update i copies j's opinion, whereas in reverse node 
update i copies i. 

By applying the methodology proposed above we find 
that, in contrast to the result for link update rules, chang- 
ing the rewiring scheme brings about different scaling 
laws for consensus time when node update rules are used. 
The results are summarized in Table HI (see supplement 
for details). With rewire-to-same, both versions of node 
update rules exhibit linear scaling of Tc , but for rewire- 
to-random one the growth is either logarithmic or close 
to exponential, as was previously observed in simulations 
in [3. 

As shown in the Table, the growth of Tc is slightly 
slower than exponential in the case of node reverse up- 
dates. We can test this prediction by considering dense 
networks in which the average degree scales with the 




FIG. 4. (Color Online) Consensus time in other model specifi- 
cations, (a) For link-update rules with moderate k, consensus 
time does not depend strongly on the rewiring rule. Theoret- 
ical predictions from Eq. (|10p fsohd line) and Eq. (|14p f dashed 
line) are compared to simulation results averaged over 1000 
networks for rewire-to-same (orange circles) and rewire-to- 
random (green squares), (b) Linear scaling for node-direct 
updates in dense networks with k = A'^/10. The theoretical 
prediction of Tc ^ N (solid line) is compared to simulation 
results averaged over 100 samples (purple circles). Parame- 
ters: 4> = 0.5, iV = 1000. 





Rewire-to-same 


Rewire-to-random 


Edge 

Node-direct 

Node-reverse 


Tc '^ N 
Tc ^ N 
Tc ^ N 


Tc -^ N 

Tc ^k log(iVfc) 



TABLE I. Summary of the scaling laws (in large A'^ and k) 
found for the mean time to consensus in various specifications 
of the adaptive voter model. 



number of nodes according to fc = cN . Here the theory 
predicts a return to linear growth, which is confirmed 
numerically in Fig. |4}d. 

We have formulated an analytical theory for the emer- 
gence of consensus in an extension of the adaptive voter 
model. By including the simple and sociologically plau- 
sible conflict resolution mechanism of seeking another 
opinion, we have shown how the formation of consensus 
may be enhanced or suppressed. This effect is manifested 
in trichotomy of scaling laws for consensus time, between 
exponential, linear and logarithmic. We also applied the 
proposed method to several other specifications of the 
model, showing how the previously observed sensitivity 
to model specification breaks down for highly connected 
networks. In all of the models investigated, consensus 
formation is driven by the intrinsic noise arising from the 
microscopic dynamical rules of the system. This noise is 
a universal feature of network models, being a result of 
the of the discrete nature of the individual interacting 
nodes and edges. We expect that the methods developed 
in this work will provide insights into emergent phenom- 
ena in other network systems. 

Acknowledgements TR thanks Alan McKane for brief 
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Supplementary material 



Derivation of the Pair Based Proxy 

In this section we describe the derivation of a low-dimensional Markov jump process which captures 
the essential features of the adaptive network model. The method involves (i) choosing a small 
number of appropriate summary statistics (for example, the numbers of certain types of edge), (ii) 
enumerating the possible events that bring about a change to these quantities, and (iii) computing 
the probability of these events and the average change resulting. 

First, we recall the model definition. Consider a network of A^ nodes and K edges in which both 
the edges and the node opinions are initially assigned uniformly at random. In each time step an 
edge («, j) is selected at random. If the edge is active (that is, the end points have different opinions) 
then with probability the edge is rewired, with node i cutting the edge and creating a new one 
to another node, selected at random, but with the same opinion as i. If the edge is not rewired 
then a third opinion is sought: another node is selected at random from the rest of the graph and 
with probability p both i and j adopt the opinion of the new node, otherwise, i and j both take the 
opposite opinion. 

The state of the system at a given time is summarized by the state vector O, whose components 
are X = [A] — [5], F = [AA\ — [BB] and Z = [AB]. There are four basic events which can alter the 
state of the system, which we will consider in turn. First, a node in state A may rewire one of its 
edges from a i? to an A. This results in [AA] increasing by one, and [AB] decreasing by one. The 
changes to the state vector 17 are captured by the vector 




s,= \ 1 I . (15) 

The probability of this event is given by 

From left to right, the contributing factors are (1) the chance of selecting an {A,B) edge from the 
K available, (2) the probability of rewiring and (3) the probability of picking the A node to rewire. 
With the same probability, a B node may also choose to rewire, resulting an increment to [BB] and 
decrement to [AB]. Similarly, the jump vector for this event is 

s,= \ -l] . (17) 




The other possible events are A or B nodes altering their states after seeking a third opinion. 
Recall that an A node will change to B with probability p if the second opinion was -B or (1 — p) 
if it was A. Taking all the factors together, the probability per time step for a node to change its 
opinion is ( 1 - 0) (p [5] + ( 1 - p) [A] ) [AB] /NK for A nodes and ( 1 - 0) (p [A] + ( 1 - p) [B] ) [AB] /NK 
for B nodes. 

To quantify the typical change to the system resulting from such an event requires some thought. 
Suppose an A node copies the opinion of one of its B neighbors; we can immediately deduce that 



[B] and [BB] will increase by one, while [A] and [AB] will decrease by one. In addition, we must 
consider the other neighbors of the node whose opinion changed. On average an A node has 



kA 



2[AA] + [AB] 



neighbors, of which we expect 2[AA]/[A] to be of type A and [Ai?]/[A] to be of type B. The node 
in question already has one neighbor that we know about, so the average number of additional A 
neighbors it has is (1 - 1/A;a)2[AA]/[A] and B neighbors is (1 - l/kA)[AB]/[A]. 

Putting this information together, we arrive at the following vector describing the typical change 
to n brought about by an A node changing its opinion: 



/ 



S3 



l-( 



-2 

-kA 

. _ J_\ 2[AA]-[AB] 
^ kA / [A] 



To write S3 in terms of X, Y and Z, we note the following: 

[A] = \{N + X), [5] = i(iV-X), 

[AA] = ]^{K + Y-Z), [AB] = Z, [BB] = ^{K -Y - Z) 

K+Y , K-Y 

kA = 2 -——- , kB = 2- 



(19) 



N + X 



N-X 



Thus 



/ 



S3 



-2 

-, K+Y 

N+X 



\-l + {K + Y-2Z){ 



N+X K+Y 



\ 



'+y)J 



(20) 



Finally, for B nodes changing their opinion, the update vector S4 is obtained in exactly the same 
way: 



S4 



V 



2 

ks 

. _ J_\ 2[BB]-IAB] 
^ kB [B] 



-1 



q K-Y 

{K~Y-2Z){j^ 



X 



K-Y ) 



(21) 



There is an implicit assumption in the above derivation. When we estimated the number of ad- 
ditional A and B neighbors belonging to the updated A node, we used the statistics of typical A 
nodes. However, we know that the node in question already has at least one B neighbor; in certain 
circumstances this fact can mean that the node is actually rather unusual (for example if there are 
not many AB edges). This assumption is essentially equivalent to the pair-approximation, which is 
regularly employed in studies of adaptive networks. The approximation is uncontrolled in the sense 
that there is no rigorous bound on the error committed, however, it is known to work well in the 
voter model provided k is not too small and is not too large 17 . 
For notational convenience, we package the update vectors into a stoichiometric matrix 5 = (si : 



S2 : S3 : S4), and collect the event probabilities in a vector 



/ 



[AB] 
K 



0/2 
0/2 



\ 



(l-0)(p[5] + (l-p)[A])/iV 

V(i-0)(p[A] + (i-p)[5])/ivy 



/ 









\ 



(l-0)(l + (l-2p)x) 
V(l-0)(l-(l-2p)x)/ 



(22) 



Here we have introduced the rescaled variables x = X/N y = Y/K and z = Z/K. For A^ ^ 1, 
application of Kurtz' theorem provides the mesoscopic equations 



d 
di 



X 



Sr 



1 



vit) 



(23) 



where 



and 



S 



/ -2 2 

1 _i -fc(^) k{^) 



with the noise correlation matrix B is given by 

k 

Slow manifold projection 



(24) 



(25) 



(26) 



As explained in the main text, we are interested in the behavior of the model when p ^ 1/2. We 
begin with an analysis of the case that p is exactly equal to 1/2. In this case the equations describing 
the system are, at last, short enough to be comprehensible: 



dx 1 
dy 



X — y 1 , , 



(27) 



In the deterministic limit N ^ 00 the system possesses a pair of slow manifolds of fixed points: the 
plane z = and the parabola 



dt 1- X^ y/N 

dz 2z ( , 2^(1-0) k{\ - <\))(\ - xy - 2z)\ 1 






(28) 



where 0^, = (A; — '2)/{k — 1). Whilst every point of the plane z = is unstable, the one-dimensional 
manifold described in f l28|) is attractive. We compute the Jacobian matrix at a point (x, y, z) chosen 



from the parabola: 



/ 



J 



\ 



I x{l+<t>-2<t>,){<t>,-<t>) x{2{<f>,~4,~k{l~<l>)){4,,- 





2-d>, 



(29) 



since it is lower-triangular, the eigenvalues of J may be read from the diagonal: 



Ai = 0, A2 
The corresponding eigenvectors are 

Vi = 

where the constants are given by 

/^i = -1 + 

fi2 = -I - 



A. 



2- 



(30) 




,V2 




.V3 




(31) 



1 



2 



(32) 



(A:-1)(A;-2)(1- 



The second two eigenvalues are negative, and large in comparison to Ai = 0, meaning that trajectories 
close to the parabola collapse quickly in the directions of V2 and v^. 

We intend to make this intuition formal by writing a reduced one- dimensional equation for motion 
along the slow manifold. We reason as follows: the full stochastic equations ( 127|1 can be thought of 
as governing the random motion of a particle x, whose coordinates tell us the state of the network 
model. Suppose that at time t the particle is sitting on the slow manifold (jS]). In a short time window, 
the noise terms in ( 1271) result in a random perturbation to x by the addition of a small 'kick' dx 
(arrow (a) in Fig. [S]); the entries of this vector are Gaussian random variables with covariance matrix 
B. Once the particle departs the slow manifold, the deterministic part of the dynamics becomes 
important again, causing a rapid decay along the direction of the fast eigenvectors (arrow (b) in 
Fig.EJ. 

We approximate the effective motion along the manifold by computing the projection oi x + dx in 
the direction of V2 and v^ onto the line through x tangent with Vi. The final position can be written 
as a; + dx', where x' is a modified Gaussian random vector with covariance matrix B' = PBP^ . 
The projection matrix P has the form P = Vi{u^Vi)u^, where u is any vector orthogonal to i>2 and 

In the present case the P has the simple form 



P 




(33) 



The zeros in the second and third columns mean that, in the neighborhood of the slow manifold, it 
is the X coordinate which dominates the dynamics. The reduced equation for x alone is 



dx 
Iti 



2(0,-0)(l-a;2) 



N 



m, 



(34) 
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x+dx 



x+dx' 




FIG. 5. Illustration of the projection. 



where ^ is a standard Gaussian white noise variable. Applying the same method when p ^ 1/2 
yields the more general equation 



^ = (l-2p)(0,-0)x(l-x2) + 



Time to consensus 



0)(l-a;2) 



N 



m- 



(35) 



In the space of {x,y,z) variables, the consensus states are the points (1,1,0) and (—1,-1,0), cor- 
responding to all nodes having opinion A or B, respectively. These are also the endpoints of the 
parabolic slow manifold (|5]). Moreover, the states x = 1 and x = —1 are absorbing boundaries in 
the reduced equation f p5|) . The question of how long it takes the network to reach consensus is thus 
mapped onto the more accessible problem of determining the time until a particle obeying the SDE 
f l35|) is absorbed. This problem is solved by textbook [l5[ methods, which we now review. 
Consider a general one-dimensional Ito SDE 



^ = f{x) + g{x)v{t). 



(36) 



We write T{x) for the expected time taken for a particle to hit either —1 or 1, given that it is at 
position X at time t = 0. It can be shown that T{x) satisfies the second-order ordinary differential 
equation 



with boundary conditions 



/W|r(.) + l,(.f^,T(.) = -1 



T(-l) =T(1) = 0. 



(37) 



(3^ 



2 p f(v) fly 

This equation is solved with the introduction of the integrating factor il){x) = e a(y)^ ■ ^ the general 
result is 



T{x) = 2 



X 1 

-1 v(y) 



dy 



1 1 rv V'(z) 

X 'ip{y) J-l g{z)'^ 



ry V'(z) , , 



1 1 

X tp{y) 



'^yj (/-I ik) I-i S$''"'V 



/-I ^^y 



(39) 
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In the special case that ip is even and x = 0, the denominator cancels with the leading factors in the 
numerator, to give 

;dzdy , 



^(y) Jo 9{zf 

where the second line comes from using the symmetry of ip about the origin to infer that the 
contribution to the integral from below the y axis is equal to the contribution from above. Applying 
this formula with / and g taken from equation (1531) yields 

nO) = j——: / -. -dzdy. (41) 

A solution to this integral may be written using hypergeometric functions, however it is not especially 
enlightening. We will instead focus on the asymptotic scaling behavior as N grows large. As discussed 
in the main text, we find three different regimes, depending on the value of e. 
In the simplest case of e = we obtain the exact solution 

For negative e we observe that if N is large, the integral over z is almost completely dominated by 
the contribution from the neighborhood of 2; = 0. The standard approach in a situation such as this 
is known as Laplace's method, or the Gaussian approximation: the integrand is well-approximated 
by a Gaussian curve centered at zero. Thus, 




(43) 



el^N' 



where the final line comes expanding the exponent —eNy"^ around its maximum at y = 1. 

For positive e, the reverse of the above logic applies: integral over z is dominated by the boundary 
z = y. Unfortunately, the contribution from the denominator diverges as ?/ — )■ 1, meaning that the 
double-integral is not accessible by Laplace's method. We change strategies at this point and instead 
focus our attention on the dynamical equation itself. 

If e < then (135|) is deterministically unstable at x = 0, meaning that after an initial perturbation 
the system approaches consensus rapidly. In this case, the role of intrinsic noise is simply to provide 
the initial perturbation. The short-time evolution of a trajectory starting at zero is well-described 
by a linear approximation. Expanding to first order around a; = we have 



| = 2.(^.-^).+ ^H^„W. (44) 
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This is an Ornsein-Uhlenbeck process, which admits an exact solution. At general times, x has a 
Gaussian distribution with variance 

{x^) = J^ (e^^(^*-^)* - 1) . (45) 

Before the deterministic part of the dynamics can take over and drive the system to consensus, we 
must wait for x to reach a magnitude of order one. Putting (a;^) = 6 and solving for t gives 

log(l + e6N) log(iV) 



4e{<p, - 0) 4£(0, - 0) + ^^^^ 
Alternative model specifications 



(46) 



In the model studied in the previous sections, we chose edges at random to update and in a rewiring 
event the new node was chosen to have the same opinion as the initial node. Both of these choices 
are somewhat arbitrary, and may be changed. In the literature, three possible update rules have 
been considered: 

Link: An edge is chosen at random. If the edge is not rewired, then one endpoint copies the other, 
chosen at random. 

Node-direct: One node is chosen at random, the other node is selected at random from the neigh- 
bors of the first. If the edge is not rewired, then the first node copies the opinion of the second. 

Node-reverse: One node is chosen at random, the other node is selected at random from the 
neighbors of the first. If the edge is not rewired, then the second node copies the opinion of the 
first. 

There are also two possible choices for rewiring scheme: 

Rewire-to-same: The new node is chosen to have the same opinion as the rewiring node. 

Rewire-to-random: The new node may have any opinion. 

The procedure developed for the LS (link update, rewire to same) specification of the model also 
applies to the other variants, as we now describe. For clarity, we are considering only the case 
p = 1/2, that is, the original adaptive voter model. 

In each of these schemes the basic events which change the state of the system remain the same, 
and thus the stoichiometric matrix S is unchanged. The rates, however, are different in each case. 
For edge dynamics, we have 

/ 6 \ 



Tls 



2K 



\ 



Tlr 



2K 



J 



f (t)[A]/N\ 
(l>[B]/N 

(1-0) 
(1-0) 



(47) 



The extra factors of [A]/A^ and [B]/N in the rewire-to-random case come from the probability of 
choosing a node of the same opinion to rewire to. For the node-direct (D) and node-reverse (R) with 
rewire-to-same, the rate vectors are 



/ 



ros 



[AB] 

N 



(p/kA \ 

(p/ks 
I - 4>)/kA 
I - 4>)/kB I 



( 



rRs 



[AB] 

N 



(p/kB \ 
(f)/kA 
1 - (p)/kB 

1 - <P)/kA J 



(4J 
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The overall denominator is N in this case since we are choosing from the population of nodes, 
rather than edges. The factors of l/k^ and l/lfcs are the result of choosing from the neighbors of 
a node; these factors are reversed in the node-reverse dynamics. Finally, for node dynamics with 
rewire-to-random we find 



roR 



N 



I <P[A]/NkA \ 
<P[B]/NkB 

(1 - <P)/kA 
V (1 - ct^)/kB j 



rRR 



N 



( cl>[A]/NkB \ 
(t)[B]/NkA 

(1 - <P)/kB 

V (1 - 0)/A;^ j 



(49) 



We begin our investigation of consensus time by examining the choice of link update rules. From 
the work of the previous sections, we know that 



T, 



LS 



N 



log(2) 



(50) 



For rewire-to-random, the SDEs for x, y and z are 

dx 1 
'dt~ 

dy_ 

dt 
dz 
Iti 



Vx{t) 



2z(l — (t))(x — y) X 1 

^ , '\ + Y + -^Vy{t) 

1-x^ k Jn 



(51) 



2z{l - <p){l - xy - 2z) 4^(1-0) 30-4 



1 



l-x2 



k 



Vz{t) 



k(l-y^) 

where the noise matrix B is given by the usual formula. There is once again a manifold of fixed 
points, although now it does not have such a simple form: 

y = x 



X 



2k{l - 1/0) ' 
(1 - 1/2) f (4 - 30)(1 - x2) - 2k{l - 0)(1 - xy) 



4(l-x2-A;(l-|/2))(i 



(52) 



From here the calculation proceeds exactly as in the case discussed in the previous sections: we 
compute the Jacobian, and project onto the slow manifold along the fast direction. The result is a 
one-dimensional approximation for motion along the manifold, namely 



dx 
'dt 



' (1 - a;2) (2fc(l - 0) + (1 - x)x(t)){x{l + x)(f> - 2fc(l - (j>)) (4 - 2fc(l 



(3-a;2), 



fc(4fc2(l 



- .t2 (1 - a;2) 02 + 4A:(1 - 0) ((1 + x^) (j) - 1)) 



■m, (53) 



where ^(t) is a Gaussian white noise variable. Exact computation of the time to absorption at 
a; = — lora; = lis difficult, however, to ffist order in large k we have 




(l-x2)0 



m 



(54) 



Applying fl39l) . we find 



Tlr 



N\og{2) ^ iV(log(2) + 0/4) ^ ^^^_2^ 



1- 



k{l- 



(55) 
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For large k, then, the choice of rewiring scheme makes very httle difference to the time taken to 
reach consensus. We test this conclusion numerically in Figure 4a of the main text. 
We move on now to study node update rules with rewire-to-same. For node-direct, we have 

dx _ 2z{l - (l)){y - x) 

— 2 VxV') 



dt 1 — y^ 

dy 2z(f){x — y) 
dt k{l — y 



2^ 



+ Vy{t) (56) 



d^_ [ Az{l - xy){l - ct>) _ l + {l + 2{k + l)z){l-ct>)-xy{2-ct>) ^ .. . ^ , ^ ^ .,^ 
dt~ \ k{l-y^Y k{l-y^) +^ <P)]+nz[t). 



and for 


node- 


reverse 








dx 
dt 


2z{l 


-(j)){x-y) 


+ r]x{t) 






dy 


2z{x 


~y){2k{l- 


■<p){l-xy)- 


{l-x- 


'^'^Kvit) 


dt 




k{l- 


-x^){l-y^) 




1 rjy[i) 


dz 


2z{l ■ 




lz-l){k-l) 


2-0/ 
1-0V 


l-xy) 


dt 




k{l~y 


'^) 










, {i-^y?- 

(1- 


-2(2^- 
^x^){l- 


-l){l-xy) 1 




- y2) 1 _ a;2 



(57) 



Although these expressions are somewhat different from each other, and from those of the rewire-to- 
same model with edge dynamics, they have exactly the same parabola of fixed points as discussed 
above, namely 

y = x, z = -[I — X ) . (58) 

In both cases, the eigenvalues on the slow manifold are 

The details of the noise projection are slightly different for each version of the dynamics. For node- 
direct we find the reduced equation 



f-f'^-')'^--^'0 ^.a-V)tV l^'"- •'"> 



and for node-reverse 



Both of these are only slight modifications of the result for edge dynamics with rewire-to-same. The 
time to consensus is found in the usual way, giving 

A^log(2) _ iVlog(2) 
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The conclusion we can draw from this exercise is that, under the rewire-to-same rule, the choice of 
edge or node dynamics has very little difference on the mean time to reach consensus. As we will 
now see, this picture is dramatically different for the rewire-to-random rule. 
For node-direct dynamics with random rewiring, the SDEs are 

dx 2z{l — (f)){y — x) 



dt 1-1/2 

dy_ ^ z(t){xi2-xy) -y) 
dt k{l-y^) 

dz z(l){2xy — 1 



X 



k{l-y^) 
and for node-reverse 

dx 2z{l — 0)(x — y) 
It ~~ 



+ Vy{t) 

-2z(l- 



(63) 



{l-xy + {k + l)z) A{z-xyz) \ 



k{l-y') 



k{l-y^f J 



Vx{t) 



1-2/2 
dy (f)zy{l — x"^) ^ 4(1 — (f)){y — x){l — xy) 
lit 



dz 
li 



k{l-y') 

(j)z{x'^ — 1) 
k{l-y') 



+ 



+ riy{t) 



+ 2z(l-0) 1 



x2)(l — y"^) 

2 2{{k + l){z - I) + xy) 



(64) 



1 — a;2 



k{l-y^) 

A{l-xy){l-z ) 
[l-x^){l-y^ 



+ 



+ Tl^{t) 



The behaviors of these versions of the model are fundamentally different from those previously 
considered: they do not posses a non-trivial manifold of fixed points. Instead, the fixed points for 
both ( l63l) and (IMll are the plane z* = 0, and the point 



•A/*: 



0, y* 







2(1-0) (2-0, 



(65) 



where 0, = {2k — 4)/(2/c — 3). The time to reach consensus in each of these systems is dominated 
by the behavior in the neighborhood of this fixed point. 



Ai 



4-30. 



A 



2,3 



2k 



2k{l 



+ (j)±J {2k{l - 0) + 0) + 8k{l - 



(66) 



Notice that for all positive k and G [0, 1] we have 

2A;(l-0) + 0< J{2k{l- 



+ 



+ 8k{l 



So, for < 0, we have that Ai and A2 are real and negative, but A3 is positive, meaning that the 
fixed point is unstable in the direction of the eigenvector 



■"3 



^ 2k{l 



{2k{l 



+ 



8k{l 
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Typically, however, this instability is weak. Expanding in large k, we find 

Ai^-2(l-0), A2^-(l-0), A3^^. (67) 

For moderate values of k then, typically trajectories collapse quickly towards (a;^,,|/^,, 2;*) in the 
direction of Vi^2; whilst perturbations in the direction of v^ are weakly amplified. We exploit this 
observation to gain analytical traction by once again projecting trajectories, this time onto the 
manifold {{x^,,y^:,z^,) + sv^ : s G M}. The resulting one-dimensional model is 

^ = \3W + bm, (68) 

where w describes the distance of the system from the fixed point, along the direction V3, b is the 
amplitude of the projection of the noise matrix in the direction V3, and C,{t) is a standard Gaussian 
white noise. In the present case we have 



Since A3 is positive the mean square deviation of w will grow with time, specifically, 

M = ^(e^^3*_i). (70) 

The time taken for \w\ to reach a magnitude of order one (i.e. the time it takes the system to break 
its symmetry and then achieve consensus) thus scales like 

For node-reverse dynamics, the eigenvalues at the fixed point are 

2(6^ - 6) 



Ai 



4-30. 
'2k 



A2,3 = -|^ ( 2k{l - 0) - ± J{2k{l - 0) - 0)' - 8A;(1 - 



(72) 



This time, 



2k{l - 0) - > J{2k{l -(f))- (py - 8k{l - 

So, for (f) < (f)^ all three eigenvalues have negative real part and the fixed point is therefore stable. To 
reach consensus requires a sufficiently large noise-drive excitation for the system to escape the pull 
of this attractive fixed point; once again a linear theory will be sufficient to describe this process. 
Expanding for large k we find 

Ai^ -2(1-0), A2^-(l-0), As^-^, (73) 

so we are again in the situation that collapse in the direction of the first two eigenvectors is very 
rapid. Projecting onto v^ as usual, we obtain the same linear stochastic equation as above: 

^ = \3W + bm- (74) 



17 
This time A3 < 0, so we may apply the same calculation as in (H3|) to deduce 



1 / AgX / 7r62 1 / iV0 \ /l67rA;3(l 



rc«-exp(-— ) ^/--^f^-exp( ^^-^^ -]^l ^^^ . (75) 



